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Abstract 

This is a first step in counting the number of multiple solutions 
in certain glassy random matrix models introduced in refs. pQ. We 
are able to do this by reducing the problem of counting the multiple 
solutions to that of a moment problem. More precisely we count 
the number of different moments when we introduce an asymmetry 
(tapping) in the random matrix model and then take it to vanish. 
It is shown here that the number of moments grows exponentially 
with respect to N the size of the matrix. As these models map onto 
models of structural glasses in the high temperature phase (liquid) this 
may have interesting implications for the supercooled liquid phase in 
these spin glass models. Further it is shown that the nature of the 
asymmetry (tapping) is crutial in finding the multiple solutions. This 
also clarifies some of the puzzles we raised in ref. j2j. 

PACS: 02.70.Ns, 61.20.Lc, 61.43.Fs 

1 Introduction 

Random matrix models can be used very effectively as simple mathematical 
toy models where many new ideas in physics, biology and economics can be 
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tested analytically ref. (UGHE]- Here we try to understand the idea of tapping 
and counting, well studied in the context of granular media, in the glassy 
random matrix model introduced in ref. pQ. There it was demonstrated that 
the matrix models with gaps in their eigenvalue distribution had multiple 
solutions and were related to the high temperature phase of certain p-spin 
glass models ref. [7j. We approach the problem in much the same spirit as 
done for spin systems in ref. jS]. This is a first step in understanding what 
happens when we tap the model ie introduce a perturbation and remove it. 
This enables us to count the number of different configurations. Studies to 
understand the fluctuation-dissipation relations and the relations between 
the dynamical and Edwards temperature in the dynamical matrix models 
awaits further work. This study will also help us understand some of the 
puzzles that we raised in ref. One of the puzzles in these models is that 
the long range correlators found in ref. JU] by mean field calculations differ 
from that found in ref. |13| |2] using the orthogonal polynomial methods. A 
resolution of this has been suggested in ref. JT] where it is claimed that the 
difference arises due to discreteness of the number of eigenvalues for double 
well models with equal depths. Here we try to understand these results using 
the method of moments. 

Most of the studies and applications of matrix models correspond to eigen- 
value distributions on a single-cut in the complex plane where the eigenvalue 
density is non-zero ref. [1] . Here we study a one hermitian matrix model with 
a more complicated eigenvalue structure. These have found applications in 
two-dimensional quantum gravity, string theory, disordered condensed mat- 
ter systems, superconductors (with complex vector potential and with impu- 
rities) and glasses. Here we study these models with applications to glasses in 
mind as discussed in refs. [T]. To illustrate some of the generic properties we 
study a one hermitian matrix model with two cuts for the eigenvalue density. 
One of the important differences observed in these models is that they have 
multiple solutions which show up in certain correlation functions. Here we 
count the number of multiple solutions and explore the possibility that these 
multiple solutions arise by taking different paths in phase space ( each path 
may correspond to a different metastable glassy state ). It is important to 
establish the correspondence between the multiple solutions and metastable 
glassy states. The barrier heights corresponding to these various solutions 
are also future goals. 

I will discuss here the matrix model with double-well potential the M 4 
model (in the Gaussian Penner model where similar things happen will be 
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pursued elsewhere). A tapping is introduced which corresponds to coupling 
the matrix model to an external source. The limit of taking the external 
sources to vanish gives different values for the moments in these models. 
This may result in different values for the partition function and hence the 
free energy. Taking different tappings corresponds to exploring the full space 
of configurations. Here we present the first steps in counting the number of 
different configurations and find it to be exponentially large. 

After this work was completed we find that in a different context results of 
exponentially large number of minima have been reported in a renormalizable 
matrix potential with Sn using a different method by Soljacic and Wilczek 
ref. 0. 

2 Notations and Conventions 

Let M be a hermitian matrix. The partition function to be considered is 
Z = J dMe- NtrV{M) where M = N x N hermitian matrix. The Haar 
measure dM = UZi dM uUi<j dM^fdM^f with M {j = Nl\f + iNl\f and 
N 2 independent variables. V(M) is a polynomial in M: V(M) = g\M + 

(g 2 /2)M 2 + (g 3 /3)M 3 + (# 4 /4)M 4 + The partition function is invariant 

under the change of variable M' = UMW where U is a unitary matrix. 
We can use this invariance and go to the diagonal basis ie D' = UMU' 

such that D' is the matrix diagonal to M with eigenvalues Ai,A 2 , A at- 

Then the partition function becomes Z = C l\f=i dXiA(X) 2 e- N ^=i y(Al) 
where A (A) = Yii<j W ~ Aj| is the Vandermonde determinant. The integra- 
tion over the group U with the appropriate measure is trivial and is just the 
constant C. By exponentiating the determinant as a 'trace log' we arrive at 
the Dyson Gas or Coulomb Gas picture. The partition function is simply 
Z = CfZoUlidXie- 3 ^ with S(X) =NJ2l 1 V(\ i ) -2V,.,,,/" A, - A,|. 

This is just a system of N particles with coordinates Aj on the real line, 
confined by a potential and repelling each other with a logarithmic repulsion. 
The spectrum or the density of eigenvalues p(x) = Y,iLi ${ x ~ K) is i n the 
large N limit or doing the saddle point analysis just the Wigner semi-circle for 
a (Gaussian probability distribution for the eigenvalues) quadratic potential. 
The physical picture is that the eigenvalues try to be at the bottom of the 
well. But it costs energy to sit on top of each other because of logarithmic 
repulsion, so they spread, p has support on a finite line segment. This 
continues to be true whether the potential is quadratic or a more general 
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Figure 1: (a). The confining potential (b). The density of eigenvalues 

polynomial and only depends on there being a single well though the shape 
of the Wigner semi-circle is correspondingly modified. For the quadratic 
potential the density is p(x) = ^\J(x — a){b — x) where [a, b] are the end of 
the cuts. See Fig. [TJ 

On changing the potential more drastically by having two humps or wells, 
the simplest example being a potential V(M) = — ^M 2 + f M 4 , the density 
can get disconnected support. The precise expressions for the density of 
eigenvalues are as follows: 



p(x) = —xJ(x 2 — a?){b 2 — x 2 ) a<x<b 

= -b<x<-a (2.1) 

where a 2 = ^\\lA — 2y/g] an d b 2 = ^[\p\ + and |//| > 2y/g, which 

is the condition that the wells are sufficiently deep. The eigenvalues sit in 
symmetric bands centered around each well. Thus p has support on two line 
segments. As \p\ approaches 2-Jg, a — > and the two bands merge at the 
origin. The density is then 

= -M<x< M 

V 9 V 9 V 9 

= otherwise. (2.2) 

The phase diagram and density of eigenvalues for the M 4 potential is 
shown in Figs. 121 
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Figure 2: (a). The double- well potential (b). Density of eigenvalues (c). 
The phase diagram 

The simplest way to determine p(z) explicitly is to use the generating 
function F(z) =< jjTrj^ > and its saddle point or Schwinger-Dyson 
equation also known in the mathematics literature as the Riemann-Hilbert 
problem F(z) = \[V'{z) + \TK{z)} with A(z) = V\zf - Ab{z) and b(z) = 
gz 2 + ji + g < jjTrM 2 > (see ref. [9J. The density p{x) is then determined 

by the formula p(z) = — ^/m^A(z). In what follows the matrix model is 
tapped ( that is a small perturbation is added which breaks the Z2 symmetry 
) and the number of solutions corresponding to the different moments of the 
model is counted. 



3 Introducing Asymmetry (tapping) 

Let us put a matrix source A, with eigenvalue a n , which will ultimately vanish 
in the partition function 

Z N (A) = J dMe -NTr(V(M)-AM)_ ^ 

Using Harish-Chandra-Itzykson-Zuber 

Zn(A) = f ft d\i^r\e- £><**)--.*) (3.4) 

where 



A(A) = Det\j-\ 



(3.5) 
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Figure 3: The Asymmetric Potential V(x) 

Then in terms of the moments the partition function becomes 

Det(m n (a k )) , . 

Z N (A) = a(q) (3.6) 

with 

m n (a) = / dxe- N[v{x) - ax] x n . (3.7) 



Let us consider m n (a) if iV goes to infinity before a — > 0. 

(I) . First take a non-z? 2 symmetric V(x) with two wells see Fig. |3] 

(i) . The saddle-point is solution of V'(x) = a see Fig. HJ 

(ii) . If a is positive we have three solutions but the action is lowest at £3. 

(iii) . £3 is still the leading saddle-point solution for a < 0. 

Therefore the behaviour of m n (a) for small a is independent of sign a. 
This corresponds to the case studied in ref. where the difference in the 
depth of the asymmetric wells is large. 

(II) . However if V is symmetric, example: V(x) = -^-x 1 + |x 4 , when 
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Figure 4: Derivative of the Asymmetric Potential V'(x) 

a — > the saddle-points are 

* c = ±-l + | + 0(a*) 

(xkO has a higher action) then 

2 # 

The integral m n is dominated by 

1 a r 
x = H 1 — for a > (J 

v 7 ^ 2 

1 a „ 

= = + - /or a < 0. 

v 7 ^ 2 

The moments are thus given by 



V(X)| 




Figure 5: The Asymmetric Potential V(x) With Two Wells Of Equal Depths 



For n even the two results are the same; but for n odd we get opposite 
signs. Note that the Z 2 symmetry would say that m n = for n odd and 

N 

a — > 0. The set of moments would be 2~ corresponding to the number of 
different possible moments (only the odd moments are different for different 



(III). We have to check whether the non-uniformity of the limits N — > oo, 
a — > may be present if V is non-symmetric but has two wells of equal depths. 

The same series of arguments follow through for the asymmetric potential 
with two-wells of equal depths as for the purely symmetric potential. Hence 

N 

there would be multiple solutions of the same multiplicity 2 "2" in the moments 
for this problem as well. This is the situation considered in ref. jTT], (though 
here only one of the 2~, the symmetric solution, as is refered to in ref. |9] 
was considered) and arrive at the same symmetric answer as ref. ^T] where 
they make the unequal wells equal (asymmetry tending to zero limit). 




(3.11) 



n). 
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4 First Steps In Counting Multiple Solutions 



Let us reformulate the problem in a slightly different way to enable counting 
and bring out some novel results in a form easily comparable to formulae in 
ref. j!2j . Consider the measure 

Z~ 1 exp(—NtrV(M) + NtrMA)d N ' 2 M (4.12) 

here V is an arbitrary polynomial, and A = diag(ao, a^-i) can be 
assumed diagonal. 

One diagonalizes M: if M = where A = diag(X Q , X^_i), the 

integral over Q is the usual Itzykson-Zuber integral on the unitary group and 
one finds: 



Pn (X ,Xi,..,\n-i) = Z-^ ^^^- exp (-Nj^ViX,)) . 

(4.13) 

Replacing powers of A in the Van der Monde with the orthogonal polyno- 
mials -Pfc(A) of the measure exp(— NV(X))dX. The partition function Z can 
then be expressed as: 



N-l N-l 



f n d\det(P k (\))expN £ (-V(X t ) + a^) 

J i=0 i=0 

v ,det ( [ dXP k {X)expN{-V{X) + aiX)) (4.14) 
A(ai) \J J 



A(a t 



Pn then becomes 

1 det(P k (X i )) iik=0 ... N _ 1 det(expNa l X j ) jj i =0 ... N - 



Pn(Xq, Ai, Atv_i) 



N\ det(JdXP k (X)ex P N(-V(X) + aiX)) k>l = ... N - 



i 



N-l 



exp -N V(Xi) . (4.15) 



i=0 



This formula has a simple structure. On introducing the functions F k (X) 
h[~^ P k (X)exp(-fV(X)) and G t (X) = exp (Ncn\ - fV(X)) we have 
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Pjv(Ao,Ai,...,Ajv_i) - detadA^CAJGKA))^...^! " (416) 

The matrix (/ cL\G , ;(A)F fe (A))j ) fc =0 ...jv-i has an inverse a^. Putting the 
three determinants together we get: 

Piv(Ao, Ai, Aat-i) = j^det{K{\i, \j))i,j=o...N-i (4.17) 

where 

N-l 

K(\,fx)= FkWauGfa). (4.18) 

k,l=0 

The kernel satisfies the property: 

[K*K\(\,p)=K(\,p). (4.19) 
Thus we obtain the determinant formulae 

(TV - nV 

p n (A , Ai, A n _i) = — det{K{\i, Aj))ij=o...n-i (4-20) 
for any n < N . The kernel has the form 

JV-l 

^(A,/i)= £ F k (X)FM (4.21) 

with = J2i a kiGi(f^)'i but F fc 7^ F fe . Thus K is not symmetric. In 

order to get further properties for K we consider the integral 



= J dX(Gi(X)F k (X)) l 



k=0...N-l 



dX^=lexp(N(-V(X) + a, A)) 

f dA 2 Ci X l exp(N(- V(X) + a, A)) 

i=0 



k 

-j^Y.Ci I dXX l exp(N(-V(X) + a, A)) 
v^fc i= o 



1 fc 



, Eam, (4.22) 
«fc i=0 
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mj are the moments. For symmetric potential V^(A) the above expression 
becomes ( using the expression for the moments found in the previous section 

) 



d\(G l (X)F k (X)), k= ^ 



.N-l 



a 



ki 



2tt 



1 JL 1 . jv 

CS^-^'^V^ °' <0 



(4.23) 



Summarizing 



/ = 



a k }, ai > 



Recall that x c = ±^= + | thus only for ±-^= > | the above result holds 
i.e. the integral eq. 1)4.23)1 has two values depending on whether ai > or 

< 0. Whereas for ±-^= < | the usual single well result as given in ref. 
is found. 

From the equation for K(X, p) i.e. eq. (J4.21)) which depend on the integral 
eq. ()4.23|) through a sum it may be possible that there are 2 N solutions for 
certain kernels this would corresponds to an exponentially large number of 
solutions depending on the path or different combinations of ai taken. For 
Pn(X , Ajv-i) and Z i.e. eq. (14.14)) and eq. ()4.15|) which are related to / 
through a determinant it is risky to consider the large N behavior of / before 
computing det[^xN]I- Counting at the level of K(X,p), pn(Xo, Ajv-i), Z 
and the free energy still remains an open one and needs a non-perturbative 
treatment (as shown in ref. 0). This will be pursued in a future work. 



11 



5 An Explicit Calculation Of The Integral 
Eq. (E23j) For The Double- Well Problem 

For the double- well matrix model the orthogonal polynomials are not known 
polynomials but we do know the form for the polynomial at large N ie when 
(N — n) fa 0(1). The polynomials are given by 



^n(A) 



1 



cos(iVC - (N - n)<j> + X + (-l) n ^)(A) + O(-) 



where /, £, 0, % and 77 are functions of A and are given by 



/(A) 
C'(A) 

cos 20(A) 

cos 277(A) 
sin 277(A) 
X(A) 



7T {b 2 - a 2 ) 

2A 2 
-Trp(X) 



sin 20(A) 



A 2 - 



(a 2 +b 2 ) 



(b 2 -a 2 ) 
2 

COS 0(A) 

A 

sin 0(A) 



(5.24) 



(5.25) 



Let us consider the eq. (|4.23J) with the above asymptotic ansatz for 0^ 
for large k then 



-N(V(X)-aX) 



'hhRe 



dX 



3 i( iVC _(iV-fc) 9 i+ x +(-l)fer ? ) -N{\\nX 2 -aX) 



/(A) 



h k Re J d\e N{ ^ lnfW+iC+ ^ N ' k N +i (~^ k ^-^ x2 + aX ~m) (5.26) 

Where 7^ is given by — (N — k) + |. In the saddle point approximation 
the exponent S(X) is to be minimized. The action S(X) 
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SW = ^« + < + - + aA + ^ k/(A) ,5.27) 

will have a first derivative which vanishes as shown below 



^-MA) +i (-l)^-^ + - + ^(A) = 0. 

(5.28) 

Where we have used the relation for £ in terms of p from eq. (|5.25|) . 
Solving for the density p(A) we get 



^ ) = i ( ^_ a) + ^__^_ m) + ( _ 1) ^ ( , 29) 

For the symmetric potential p(A) = -VA 2 — 6 2 , in the large N limit ne- 
glecting the last terms, as the equation for A is quadratic there are two 
solutions to the equation as shown below 



, 9 uAa a 2 — b 2 

A - -r4 7T + 



i/i 2 + l) (i/i 2 + l) 



pa ± 1^ 



^ ^ 4(a2 " 62) (5.30) 



2(^ 2 + i) 2\ii> 2 + iy ~ (i^ 2 + l)' 

Thus in the saddle point approximation the integral I for large k becomes 



J ± = J P fe (A ± )e- iV(y(A±) - aA±) + h.o. 

(5.31) 

Where Jo is a constant. Hence we have shown in an explicit example for 
the symmetric double well potential that the integral eq. ([4.23)1 for large k in 
the saddle point approximation has two solutions, which solution is choosen 
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depends on whether a > or < 0. This result indicates the possibility that the 
kernel, partition function, free energy can have 2 N solutions depending on 
the path {ai} taken as these functions all depend on the integral /, eq. ()4.23|) . 
Thus here evidence is presented that there exists an exponentially large num- 
ber of solutions, i.e. e Nln2 , in the double well matrix models depending upon 
the path taken in parameter space {a/}. It will be interesting to explore the 
possibility that these exponentially large number of solutions correspond to 
the metastable solutions of the supercooled p-spin glass that these random 
matrix models map into. 

6 Conclusions 

We have been able to map the problem of counting the number of multiple 
solutions found in ref. [Hj to a moment problem. The multiple solutions 
were discovered in the recurrence coefficients of the orthogonal polynomials 
in ref. (§]. It was known that there are an infinite number of solutions. The 
counting problem is mapped onto counting the number of ways to get differ- 

N 

ent moments. The set of moments grows exponentially as 2~ . In order to 
show this we have to introduce a small perturbation which breaks Z2 symme- 
try into the moment integral and then take the small asymmetry parameter 
to zero ( which we call tapping the matrix ). As an added bonus we are 
able to understand some of the puzzles and controversies that are found in 
ref. jU] and studied in ref. |T3J |2j . The counting at the level of the kernel, 
Ptv(Ao, Atv_i), Z and the free energy still remains an open one and needs 
a non-perturbative treatment. This will be pursued in a future work. 

The number of moments in these random matrix models are exponentially 
rising with N. These matrix models are connected with the high temperature 
phase of structural glasses as has been discussed in refs. [UEj. There could be 
interesting properties of the supercooled liquid phase which may be explored 
analytically in these simple models. For example it will be worthwhile to 
study how the metastable states of the liquid are related to the different 
paths of taking the small perturbation parameter, as introduced here, to 
zero. Future work on finding barrier heights is underway. 
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